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Abstract 

We introduce antiferromagnetic quantum fluctuations into quantum annealing in addition to 



o 

(N 

^r<. the conventional transverse-field term. We apply this method to the infinite-range ferromagnetic 



p-spin model, for which the conventional quantum annealing has been shown to have difficulties to 
find the ground state efficiently due to a first-order transition. We study the phase diagram of this 



i-C . system both analytically and numerically. Using the static approximation, we find that there exists 

P-c 



a quantum path to reach the final ground state from the trivial initial state that avoids first-order 
transitions for intermediate values of p. We also study numerically the energy gap between the 
ground state and the first excited state and find evidence for intermediate values of p that the time 



CN . complexity scales polynomially with the system size at a second-order transition point along the 

> 

^^ I quantum path that avoids first-order transitions. These results suggest that quantum annealing 



would be able to solve this problem with intermediate values of p efficiently in contrast to the case 



cn . with only simple transverse-field fluctuations. 
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I. INTRODUCTION 

To find efficient algorithms for combinatorial optimization problems is one of the impor- 
tant goals of computer science. An algorithm is efficient if its running time is bounded by a 
polynomial in the problem size. If an efficient algorithm for a problem is known, the problem 
is considered as easy. However, most of interesting combinatorial optimization problems are 
hard, i.e., even the best known algorithms take running times growing exponentially as the 
problem size increases |lj, l2| . 

Many combinatorial optimization problems can be translated into physics problems of 
finding the ground state (optimal solution) of an Ising spin system ^]^. The cost function 
corresponds to the energy of the system. This transformation enables us to study combi- 
natorial optimization problems by ideas and methods developed in physics. An interesting 
example is quantum annealing. 
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Quantum annealing [5|-l7| (and its cousin, quantum adiabatic computation [8|) is a quan- 
tum algorithm to obtain an approximate solution for a combinatorial optimization problem, 
which often outperforms simulated annealing [9|, |lO|, another heuristic algorithm coming 
form physics. Unlike simulated annealing, QA uses tunneling effect caused by quantum fluc- 
tuations to search for the optimal solution. By controlling the strength of the fluctuations 
properly, we can reach the solution with a high probability. 

To be more explicit, let us consider the following time-dependent Hamiltonian: 

H{t) = s{t)Ho+{l-s{t))V, (1) 

where Hq is the target Hamiltonian, whose ground state is an optimal solution, represented 
in terms of the z components of the Pauli matrices a^ {i = 1, . . . , N). The symbol A^ denotes 
the number of spins. The operator V is arbitrary as long as it does not commute with Ho and 
has a unique trivial ground state. This noncommutativity introduces quantum fluctuations 
into the system, causing state transitions. It is thus called the driver Hamiltonian. A typical 
example of the driver Hamiltonian is the transverse-field operator Vxp = — Yli=i ^fi where 
the (jf (i = 1, . . . , A^) are the x components of the Pauli matrix. The control parameter s{t) 
starts at zero (s(0) = 0) and increases monotonically to unity. We assume that s(r) = 1, 
i.e., the running time of QA is r. For simplicity, the linear function (s(t) = t/r) is adopted 
in most studies. We then calculate the time evolution starting from the trivial ground state 



|\I/(0)) by solving the Schrodinger equation, 

ij^mt)) = Himit)),o<t<T. (2) 

Note that H{0) = V and H{t) = Hq, and the initial state |\I/(0)) is the ground state of 
H{0). If we change the control parameter slowly (r ^ 1), the state would stay very close to 
the instantaneous ground state during the time evolution. We can then achieve the ground 
state of Hq at t = r, which is the optimal solution. 

The condition to stay close to the ground state is expressed as r ^ ^nfin according to the 
adiabatic theorem [1C|, where Amin is the minimum energy gap from the ground state. Thus 
the minimum gap determines the efficiency of QA for a given problem. In the case that the 
minimum gap decays exponentially with the system size as Amm oc exp(— aA^), the running 
time increases exponentially. This means that QA cannot solve the problem efficiently. 

Jorg et al. have shown that QA with the conventional transverse-field operator costs 
exponentially long time to reach the ground state of the ferromagnetic p-spin model for 



p > 2 ll|. The (target) Hamiltonian of this model is given by 
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In the p — )■ oo limit, this model reduces to the Grover problem [ll|, ll2|, which any known 
algorithms, even quantum algorithms, cannot solve efficiently. Jorg et al. have also shown 
that this system undergoes a quantum first-order phase transition during the time evolution, 
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which is a characteristic feature of hard optimization problems |13l-ll5|. 

Nevertheless, the above result does not necessarily suggest a complete failure of QA for 
this simple problem (|3]). Note that the above result has been derived using the transverse- 
field operator as a driver Hamiltonian. This implies that different operators may lead to 
improved performance. We show in this paper that an operator, which induces antifer- 
romagnetic fiuctuations, significantly improves the efficiency of QA for this model with 
intermediate values of p. 

We organize the present paper as follows: In Sec. HIl we define a quantum driver operator 
Vaff and the total Hamiltonian H{t). We then explain the idea of QA using antiferromag- 
netic fiuctuations. In Sec. Illlt we calculate the partition function with the static approxi- 
mation and derive self-consistent equations. From these equations, we analyze numerically 



phase diagrams for finite p in Sec. IIV[ Numerical calculations show that first-order transi- 
tions can be avoided if we change the control parameters ingeniously. In Sec. |Vl we discuss 
the p ^ oo limit. Finally, Sec. |VT]is devoted to the conclusion. 



II. ANTIFERROMAGNETIC FLUCTUATIONS 

The main proposition of this paper is that an introduction of the following antiferromag- 
netic interaction 

Vaff^+N^-Y.^:) (4) 

4=1 

in addition to the conventional transverse-field term Vtf greatly facilitates the process of 
quantum annealing as shown in the following sections. The total Hamiltonian is therefore 

H{s, A) = s{XHo + (1 - A)yAFF} + (1 - s)\/tf, (5) 

where the parameters s and A should be changed appropriately as functions of time. The 
initial Hamiltonian has s = and A arbitrary, and the final one has s = A = 1. Intermediate 
values (s. A) should be chosen according to the prescription given below. 

This idea somewhat resembles that of quantum adiabatic algorithm with different 
paths [16|. This latter approach also considers a total Hamiltonian which consists of three 
parts: a target, a driver, and another Hamiltonian He corresponding to Vaff- Whereas 
Vaff is defined uniquely, the components of He are chosen randomly according to some 
prescription. For this system, one calculates repeatedly the time evolution by changing 
He- Then, in some instances, one reaches the optimal solution. In contrast, our approach 
involves no stochastic processes and therefore only a single run achieves the goal. 



III. ANALYSIS WITH THE STATIC APPROXIMATION 

We now confine ourselves to the ferromagnetic p-spin model with odd p of Eq. (jS]) as 
the target Hamiltonian and analyze the properties of H{s,X). This section is devoted to 
analytic computations. 



A. Partition function 



We first calculate the partition function of the system of Eq. ([5]) at finite temperatures. 
The partition function can be written in the following form using the Suzuki- Trotter for- 
mula J17l |: 

Z = lini Zm 
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where "^r^zy denotes the summation over all spin configurations in the z basis, and Kcrl"}) = 
<S>i=i Wi)- The state \a^) is the eigenstate of a^, having the eigenvalue a^ (= ±1). Similar 
notations will be used for the x basis. 

We then introduce the following M closure relations: 

!(«)= E IK(«)})({^r(«)}l 

X E IK(«)})(K(«)}I, (7) 

where a = 1, . . . , M. Inserting 1(a) just before the ath exponential operator involving af 
in Eq. ([6]), we have 

M „ , ,. . AT 

P 



KH}K{°)}"=i 



/3gAiV /l ^ 



(T, [a. 



1=1 



M livi-^^^"V +^i^E^^ 

4=1 



a 



i=l 



AT 



xU(^!i<y)WK<y)){^K<y)WK<y + i)) 



(8) 



i=l 



with periodic boundary conditions such that erf (1) = af (M + 1) for i = 1, . . . , A^. 

To simplify the spin product terms (X]j=i ^iX(^)I^Y ^-^d (X]i=i ^^i.(^)/^Yi we introduce 
the following integral representation of the delta function: 
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Here, m denotes the magnetization (order parameter), and its conjugate variable is m. Using 
Eq. IQ, we can rewrite Zm as 

M 
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X J]exp[m^(a)a,^(a) +m"(a)af(a)]((T,^(a)|af(a))(af(a)|a,"(a + 1)). (10) 

We have neglected a few irrelevant constants. Since the spin product terms have disappeared, 
we can perform the summation over all spin configurations independently at each site. Then, 
we obtain 

Zm= I 'I \\dm''{a)dm''{a)drrf{a)dfrf{ 
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where the trace means the summation over the spin variables cr^(Q;), cr^(a) (a = 1, . . . , M). 
Note that the exponent in Eq. (ITT]) is proportional to N . Thus, the integrals over the 
variables are evaluated by the saddle point method, which is to take the maximum value of 
the integrand as the result of integral (see, e.g.. Appendix A.l of [l8|). The saddle point 
conditions for m^ia] and m^ia) lead to 
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m^{a) = j^{{l-s)-2s(l-\)m^{a)}. 



(12) 
(13) 



We now use the static approximation, which removes all the a dependence of the param- 
eters. We will check the validity of this approximation in Sec. IIVI After this approximation, 
we can easily take trace in Eq. (ITTl) by the converse operation of the Trotter decomposition. 



Then, using Eqs. (fT2|) and (fT3!) . we finally obtain 

dm' dm'' exp[-A^/3/(/3,s,A;m^m^)], (14) 

where /(/3, s, A; m^, m^) is the pseudo free energy defined as follows: 

/(/3, s, A; m^ m^') = (p - l)sA(m^)P - s(l - A)(m^)2 



-ln2cosh/3y {psA(m^)p-i} + {l - s - 2s(l - A)m^} . (15) 



The saddle point equations are thus 

psX{m,')'^^^ 
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^J {ps\{m^)P-^Y + {l-s- 2s(l - A)m^}^ 



X tanh/3y {psA(m^)P-i} + {l - s - 2s(l - A)m^} , (16) 

l-s - 2s(l - A)m^ 
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{psA(m^)p-i} + {1 - s - 2s(l - A)m^} 
X tanh/3A/{psA(m^)P-i}^ + {l - s - 2s(l - A)m^-}1 (17) 



B. Low temperature limit 

We next derive self-consistent equations in the low temperature limit to examine quantum 
phase transitions. Since the start of the QA process belongs to the paramagnetic phase and 
the goal is the ferromagnetic phase, a quantum phase transition inevitably occurs in the 
course of time evolution. 

It is useful to consider two possibilities separately depending on whether the argument 
of the square root in Eqs. (TT6|) and (TT7I) vanishes or not. We start our discussion from the 
latter case. 

When the square root in Eqs. (TT6|) and (TT71) assumes a finite value, the hyperbolic tangent 
tends to unity in the /3 — )■ oo limit. Then, we have 

m' = psXjmr-^ ^^g^ 



{psX{m')p-^} + {l-s- 2s(l - X)m^y 
l-s- 2s(l - A)m^ 

{psX{m^)p-^Y + {l - s - 2s(l - X)m^y 



m^ = '-'- ^''^ - "''"' (19) 



The pseudo free energy (fT5l) becomes 

/(s,A;m^m^■) = {p - l)sX{m'y - s{l - X){m'')^ 
- ^J {psX{m^)P~^Y + {1 - s - 2s(l - X)m^Y. (20) 

Equations f lTS]) and fITI?]) have a ferromagnetic (F) solution with m^ > and a quantum 
paramagnetic (QP) solution satisfying m^ = and m^ ^ 0. Substitution of m^ = into 
Eq. ^ yields 

|1 -s-2s(l- A)m^|' ^ ^ 

i.e., m^ can be ±1. However, m^ = —1 is not a proper solution since, with m^ = —1, 

l-s-2s(l-A)m^ = l-s + 2s(l-A) >OforO<s<l,0<A< 1 , which leads to m^' = 1 

according to Eq. f l2ip . The other possibility m^ = 1 satisfies Eq. (12 ip when s < 1/(3 — 2A). 

Therefore the QP phase can exist in the region < s < 1/(3 — 2A), and its free energy is 

/qp(s,A) = -sA + 2s-1, (22) 

which is independent of p. 

The free energy of the F phase cannot be obtained analytically for general p. However, we 
can evaluate it in the p — )■ oo limit as follows: In this limit, Eq. f llSp reads m^ = or m^ = 1. 
The latter solution corresponds to the F phase. The magnetization in the x direction is zero 
since Eqs. f lTSj) and f lT9|) satisfy (rra^)^ + (m^)^ = 1. Substituting the values of magnetization 
into Eq. (!20|) and taking the limit p — )■ oo, we find 

/f(s,A)|p^oo = -sA. (23) 

Let us next consider the case, where the argument of square root in Eqs. ( !T6|) and f lT7|) 
vanishes. We then assume that m^ and m^ tend to the following values as /9 — i- oo: 

such that the argument of hyperbolic tangent approaches a finite constant: 

/3^{psX{m')p-^Y + {l-s- 2s(l - A)m^}^ -^ c. (25) 

In order to find a non-trivial solution, it is also necessary to assume the following relation: 

psX{m^y~^ 



1 - s - 2s(l - A)m^ 



-^ 0. (26) 



Under these assumptions, Eqs. (fT6|) and (fT7|) read m^ = and m'^ = tanh c. These equations 
satisfy the condition flM|) if we choose c such that tanhc = (1 — s)/2s(l — A). 

Unless s = 1, the magnetizations JHM satisfy the condition of QP solution, we then call 
this phase QP2 in order to distinguish from the QP phase described before. The free energy 
of QP2 phase is obtained in the limit (^^ and /3 — > oo under the assumption f l25|) : 

The domain of applicability of the free energy (1271) is restricted by 1/(3 — 2A) < s < 1 since 
1(1 — s)/2s{l — A)| = |tanhc| < 1 and s ^ 1. This region of s will be called the QP2 domain 
hereafter. 

C. Phase transition on the line A = 

Although we cannot solve explicitly the self-consistent equations for general values of the 
parameters, it is possible to solve them in the case of A = 0. In the following discussion, we 
show that a second-order phase transition occurs at the point (s. A) = (1/3, 0) for any value 
of p. This result is independent of the target Hamiltonian. 

We start from the analysis of the phase diagram on this line. In this case, Eq. (llSp 
reduces to m^ = 0. It thus suffices to consider the QP and QP2 phases. These phases do 
not have a common domain of definition: The QP and QP2 phases are defined on the region 
< s < 1/3 and 1/3 < s < 1, respectively. Thus, the former (latter) range is the QP (QP2) 
phase, and a phase transition exists at s = 1/3. 

The magnetization in the x direction of the QP2 phase (p4|) is identical to that of the QP 
phase at the transition point. This means that the transition is of second order. The free 



energy /qp(s,0), of course, smoothly connects to /qp2(s,0) at the point |19 |. 

Though the existence of a second-order transition does not hamper the efficiency of 
QA [20], paths of QA must not follow this line because quantum fiuctuations completely 
disappear on this line: The total Hamiltonian H{s, 0) is diagonalized in the x basis. Quan- 
tum state transitions do not occur, and the system does not perform quantum annealing 
processes. 

In the following section, we will show that paths exist to avoid first-order transition in 
the region A > at least for 5 < p < 21, and possibly for any large but finite p. 

9 



IV. NUMERICAL RESULTS 



A. Phase diagram 



Let us next analyze numerically the phase diagram on the s-X plane for finite values of 
p. We construct the phase diagram as follows. We first solve numerically the self-consistent 
equations (TT8l) and (TT9|) for a given value of p and at a point (s, A) in the phase diagram and 
then evaluate the corresponding free energy. By comparing all possible solutions and their 
free energies including /qp2, we identify the stable solution having the smallest value of the 
free energy. 

It is useful to show the dependence of the free energy on s for some values of p and A as 
in Fig. [U We have confirmed numerically that the free energy lies below /qp2 in the QP2 
domain, and the QP2 phase is completely suppressed by the other phases. This system thus 
undergoes a quantum phase transition from the QP phase for small s to the F phase for 
large s. 

To determine whether the transition is first order or second order, we show the behavior 
of the magnetization m^ in Fig. |2J The parameters of the figure correspond to those in Fig. [TJ 
When A = 0.1, the magnetization m^ for p = 3 has a small jump at s = 0.3544(1), and m^ 
for p > 5 decreases continuously from unity to our numerical precision. This means that 
m^ for p > 5 increases continuously from zero to a finite value. Therefore a second-order 
transition occurs for p > 5 at A = 0.1. The same is true for A = 0.3 in the sense that there 
exists a second-order transition at the boundary of the QP2 phase for p > 5. 

A remarkable fact is that the magnetization for some parameters (e.g., A = 0.3, p = 11) 
in Fig. |2] jumps within the F phase. This discontinuity results in an exponential decrease 
of the energy gap as A^ increases. There exists a first-order transition within the F phase. 
However, this unusual behavior disappears for smaller values of A for any finite p, excluding 
p = 3, that we checked. Thus for smaller A, only a second-order transition takes place as we 
increase s from zero to a value close to unity. 

The resulting phase diagrams are shown in Fig. |3] for p = 3,5, and 11. We see that 
a boundary of second-order transition exists for small A and p > 5. It is observed that 
one can reach the F phase from the QP phase by choosing a path that avoids a first-order 
transition as long as the first-order F-F boundary does not reach the A = axis, which 
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FIG. 1. Free energy vs. s for some values of p. The parameter A is 0.1 (top), and 0.3 (bottom). 
The dash-dotted hne in hght green represents the free energy of the QP phase, Eq. ([22]), the thin 
sohd hne in blue is for the F phase, Eq. (j23p . and the thick solid line in red for the QP2 phase, 
Eq. (j27p . The vertical dashed line denotes the lower limit of the QP2 domain (s = 1/(3 — 2A)). 
Although it is difficult to discern in the present scale, all the data for finite p we studied have lower 
values than that of /qp2 in the QP2 domain. 



happens probably only in the limit p — ?■ oo as we shall see below. 
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FIG. 2. Magnetization m^ corresponding to Fig. [H The solid line represents the x component 
of magnetization of the QP2 phase ()24p . and the vertical dashed lines are the same as those in 
Fig. [TJ For A = 0.1 (top panel) and p > 5, a second-order transition occurs at the boundary of 
the QP2 domain; The magnetization decreases continuously from unity to zero. In contrast, the 
magnetization for A = 0.3 (bottom panel) has a jump. 



B. Energy gap 



We next study the behavior of the energy gap across the phase transitions found in the 
previous section. 



To calculate the energy gap for large N, we adopt the method used in ll|. The Hamil- 
tonian under consideration is expressed by the components of total spin operator S^'^, thus 

12 





1 

0.8 
0.6 
0.4 
0.2 






1 1 


F 




:_^- 


,-'■ 


- 


1 1 


QP 

1 1 



0.2 0.4 0.6 0.8 1 



0.2 0.4 0.6 0.8 1 



0.2 0.4 0.6 0.8 1 



FIG. 3. Phase diagrams on the s-X plane for p = 3 (left), p = 5 (middle), and p = 11 (right). The 
dash-dotted line represents the boundary of the QP2 domain (s = 1/(3 — 2A)), where a transition 
takes place between the QP and F phases. For large A, the QP and F phases are separated by the 
horizontal phase boundary (QP-F boundary). The thick solid line in red represents the first-order 
transition, and the thin solid line in light green is for the second-order transition. For p = 5 and 
11, the magnetization jumps on the dashed line in blue (F-F boundary) within the F phase |2l| . 



commuting with the total spin S. Since the total angular momentum is conserved during 
the time evolution, we have to pay attention only to the subspace that has the maximum 
angular momentum S = N/2. The dimension of this subspace is A^ + 1, which greatly 
enhances the possible system size to A^ ~ 100. 

It is useful to first verify the validity of the static approximation. Figure H] shows a 
representative energy gap with a second-order phase transition: As one sees in the enlarged 
view shown in the bottom panel, the gap shows wiggly behavior for a finite range. The 
wiggly behavior starts at s ~ 0.4184 for A = 0.3, which corresponds to the left end of the 
QP2 domain and also to the second-order transition point between the QP and F phases. 
The same behavior terminates at s ~ 0.4676 for A = 0.3, corresponding to the first-order F-F 
boundary. These two transition points evaluated analytically using the static approximation, 
Eqs. ( !T8|) and ( !T9|) . are shown in dashed vertical lines in Fig. H] and agree fairly satisfactorily 
with the numerical results, as A^ increases, for the interval where the gap is very small. 

The rightmost local minimum of the energy gap in Fig. H] behaves differently from other 
local minima and decays exponentially as A^ increases as shown in Fig. El This is expected 
from the jump in the magnetization shown in Fig. [2] because a jump implies a first-order 
transition though the system is ferromagnetic in both sides of the transition point. Although 
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FIG. 4. Top panel: Energy gap vs. s for p = 11 and A = 0.3. The vertical dashed lines represent 
the boundary of the QP2 domain at s ~ 0.4167 and the F-F boundary at s ~ 0.4701. The bottom 
panel is the enlarged view of the top panel for A^ = 140. 

this is not the global minimum, it will affect the efficiency of QA for much larger systems 
where the rightmost one will become the global minimum since the other local minima decay 
only polynomially as shown below. 

Figure E] shows the size dependence of local minima of the energy gap for p = 5 and 
A = 0.1. All minima shown here decay polynomially. In Fig. [7] the global minimum of 
energy gap for selected p is depicted as a function of A^ at A = 0.1. For any value of p, the 
gap closes polynomially at least up to the system size we studied, A^ = 160. 
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FIG. 5. The rightmost local minimuin of the energy gap as a function of N for p = 11 and A = 0.3 
on a semi-log scale. The gap closes exponentially with A^. 
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FIG. 6. Energy gap vs. N at local minima for p = 5,X = 0.1 on a log-log scale. We number the 
minima from left to right. No gaps vanish exponentially up to the size studied here. 

The above results suggest that first-order transitions will be able to be avoided if we 
choose a path around A = 0.1 when we reach the F phase from the QP phase by increasing 
s as long as p is not too small and not too large, 5 < p < 21. It is then interesting to see 
what happens in the limit of large p. 



15 







1.0x10 



^ 1.0x10"^ 
^ 1.0x10"^ 



LU 



1.0x10 



■12 





1 




1 


t 


« ^ 


+ 




■ 




- p=5 


+ 




'^^^^^^m^' 


: P-9 


X 




^^**^^^ 


. p=13 


X 




-^1,^^ 


-p=17 


□ 




^% 


■p=21 


■ 
1 




1 



10 



N 



100 



FIG. 7. Minimum gap vs. A^ for some values oi p on a fog-log scale for A = 0.1. Each data scales 
polynomially. 



V. PHASE DIAGRAM IN THE INFINITE-^ LIMIT 



The ferromagnetic p-spin model reduces to the Grover problem in the p — )■ oo limit UJ, 
121]. The goal of the Grover problem is to find the desired item in an unsorted database 
containing 2^ items. Whereas classical algorithms require a time of 0(2^) to find the 
desired item, the quantum algorithm, called the Grover algorithm, costs only a time of 
0(2^/^), quadratic speed-up though the time complexity still scales exponentially. 

Farhi et al. have proposed a QA version of the Grover algorithm [22|], which adopts the 
transverse field as a driver Hamiltonian. Unfortunately, the time complexity is the same 
as that of classical algorithms. However, Roland and Cerf have improved the efficiency of 
QA by adjusting the evolution rate s(t), then reproducing the quadratic speed-up, and they 
have proved that their algorithm is optimal 23|. This result indicates that our approach 
cannot avoid jumps of magnetization in the p — ;■ oo limit. It is therefore interesting to study 
how this difficulty appears in our method. 

To this end, it is instructive to study the behavior of the free energy and magnetization 
for large but finite values of p. The free energy in Fig. [1] is seen to approach the asymptotic 
values in Eqs. (1231) and (1271) from below. Hence, the QP2 phase does not appear for any 
finite p. From Fig. 12], we observe that the magnetization in the x direction is close to the 
QP2 phase magnetization (JMI), shown in red solid lines, in the region where the free energy 
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approaches /qp2- The magnetization in the z direction is 




since the QP2 phase does not appear. 

We extrapolate these results to the case of p — )■ 00. That is, while the free energies 
are described by Eqs. ( l22l) . ( 123|) . and ( 1271) . the magnetization in the QP2 phase is given by 
Eq. ( I28l) . To be precise, this is not the QP2 phase since the magnetization in the z direction 
is nonzero. With a caution on the domain of QP and QP2 in mind, we compare the values 
of the free energy of the three phases and obtain the phase diagram as in Fig. [HI The F 
phase and the QP phase are separated by a horizontal phase boundary. The boundary of 
second-order transition is given by s = 1/(3 — 2A) (A < 1/2), and the first-order F-QP 
transition boundary is s = 1/2 (A > 1/2). Solving /rlp-i^oo = /qp2, we get the F-F boundary 
as 



1-2^^^ (29) 



(2A-1)2 ■ 

The figure shows that an abrupt change of magnetization, a first-order transition, is in- 
evitable in the limit p — )■ 00. 

VI. CONCLUSION 

In the present paper, we have introduced a new approach to QA using antiferromag- 
netic quantum fluctuations. This approach adopts two types of the driver Hamiltonian, the 
transverse-field term Vtf and the transverse antiferromagnetic two-body interaction term 
V^AFF in Eq. (H. 

We have applied this method to the ferromagnetic p-spin model, which was considered 
to be hard to find the ground state for p > 2 with a simple QA [ll|. We have evaluated 
the efficiency from the phase diagram and the minimum values of the energy gap. Numeri- 
cal calculations have shown that the phase boundary of second-order transition appears for 
p > 5. However, the magnetization in the F phase jumps for large p. Although the boundary 
at which the magnetization in the F phase jumps extends as p increases, we have confirmed 
numerically that there remains a region where the magnetization changes continuously at 
least for 5 < p < 21 when A = 0.1. This indicates that QA can solve the problem efficiently 
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FIG. 8. Phase diagram in the hmit p — t- oo. Three hnes represent the same phase boundary as those 
in Fig. [3l The QP phase has the magnetization m^ = 0. The F phase above the F-F boundary, 
shown dashed in blue, has the magnetization m^ = 1 and the phase below the F-F boundary has 
< m^ < 1. 

in this case. In fact, we have calculated the minimum gap up to A^ = 160 and have con- 
firmed that it vanishes polynomially on the second-order phase boundary. Thus, QA with 
antiferromagnetic fluctuations is an efficient algorithm at least for 5 < p < 21. We expect 
to be able to avoid the difficulty of exponential complexity for larger value of p as long as it 
is finite. It is an interesting problem to study if the present method improves the efficiency 
of QA for other systems. 
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